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We demonstrate a magnetometry technique using nitrogen-vacancy centres in diamond which 
makes use of coherent two-photon transitions. We hnd that the sensitivity to magnetic fields can 
be significantly improved in isotopically purihed diamond. Furthermore, the long-term stability of 
magnetic held measurements is signihcantly enhanced, thereby reducing the minimum detectable 
held variations for both quasi-static and periodic helds. The method is useful both for sensing 
applications and as a spin qubit manipulation technique. 


Quantum interference in spin-1 systems has played a 
central role for precision measurements and tests of quan¬ 
tum mechanics using ensembles of atoms and molecules 
mn]. Coherent quantum control of three-level systems 
has furthermore been central to schemes for the coher¬ 
ent generation of single photons and the storage of light 
W- Solid-state defects have recently emerged as power¬ 
ful systems for applications in sensing, quantum informa¬ 
tion and quantum communication due to the discovery 
of emitters with long spin-coherence times in a multi¬ 
tude of crystalline hosts [3 ig. Single defects are par¬ 
ticularly promising for nanoscale sensing applications as 
they promise unrivalled sensitivity-to-size ratios. The ni¬ 
trogen vacancy (NV) defect in diamond, as well as several 
defects in silicon carbide, have shown particular promise 
as carriers of quantum information or as a sensors in bulk 
and nanocrystalline systems [THIS]. 

The electronic spin-1 ground state of the NV centre 
can be subdivided into the qubit subspaces {|0), | ± 1)} 
and {| — 1),| + 1)}. We demonstrate the characterisa¬ 
tion and manipulation of a qubit in the {| + 1), | — 1)}- 
subspace and make use of coherent interference effects to 
reduce the spin-1 system to an effective two-level system 
[n HZ]. We show that, under application of a mod¬ 
est magnetic bias field, this qubit displays significantly 
improved coherence times compared to the traditionally 
used {|0), I ± l)}-systems. Our measurements indicate 
that the reduction in dephasing rates is brought about by 
a reduced sensitivity to strain and transverse magnetic 
field fluctuations. These improvements in turn lead to 
significantly higher sensitivity for measurements of both 
arbitrary and periodic magnetic fields. The manipulation 
methods used herein furthermore enable efficient control 
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FIG. 1. a) The NV centre is driven with microwave fields 
Dm tv brought to the sample using an external microwire. A 
wire perpendicular to the axis of the NV centre carries cur¬ 
rents I used to create the arbitrary and periodic magnetic 
fields, b) Level scheme of the spin-1 ground state, shown 
with the main detuning frequencies A+,_ and driving fields 
used in the measurements. 


over the NV qutrit manifold, which have the potential to 
significantly reduce the physical requirements for large- 
scale quantum information systems [18]. We anticipate 
that our findings will be beneficial for the application of 
NV centres as sensors and in quantum technology. 

For our measurements we use a single-crystal artificial 
diamond created by chemical vapour deposition. An epi¬ 
taxially grown layer of 50 jam thickness, isotopically puri¬ 
fied to 99.999% ^^G, hosts the single NV centres [19]. The 
experimental apparatus used for the measurements is a 
standard home-built confocal microscope setup. Optical 
excitation is provided by an intensity-stabilised diode- 
pumped solid-state laser at 532 nm. We monitor the flu¬ 
orescence of the NV centres using an oil-immersion mi¬ 
croscope objective with a numerical aperture of 1.42, and 
photons are collected through a dichroic mirror using two 
single-photon counting modules. Microwave fields were 
applied through a 100 /am wire gold spanned over the di¬ 
amond surface [Fig. [^a)]. 

We manipulate the NV centre spin state using four dis¬ 
tinct microwave fields, two resonant with the |0) | i 1) 

transitions, and two fields which are detuned relative to 
the transition resonances by [Fig. [^b)]. The am¬ 

plitudes of these four fields give full control over the state 
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of the spin-1 system. The Hamiltonian of the system is 
given by 


H = 


0 11 + 0 
11 + 2 ^+ 11 _ 

0 11_ 2(A+-A_) 


(1) 


in a frame rotating at the frequency Urot = A+ — A_ 
m- 11+are the Rabi frequencies of the driving fields 
detuned by A+^_. These fields can create an effective 
coupling between the states | + 1) ^ I ~ 1): which we 
can use to characterise the dephasing properties of this 
system. We apply a bias magnetic field of 30 mT and 
can therefore neglect small transverse magnetic fields and 
strain components [5]. 

We first measured the spin coherence lifetime in the 
{| + 1), I — 1)} qubit system. For this purpose, the NV 
centre is initialized in the state |0), after which a reso¬ 
nant TT-pulse transfers the population into | — 1). From 
here, a 7r/2-pulse performed with both 11+ and 11_ brings 
the system into the superposition (| + 1) + | — l))/^/2. 
The driving amplitudes are both controlled with Gaus¬ 
sian envelopes, and the 11+ pulse is delayed with respect 
to the H- pulse. This sequence is a stimulated Raman 
adiabatic passage (STIRAP), as used for atomic samples 
in the optical domain m- We interrupt the sequence 
halfway through, at the point where the desired super¬ 
position is achieved. After a time delay the electron is 
transferred further into | + 1) using the second half of 
the STIRAP sequence. Subsequent measurements of the 
population in the three eigenstates of the undriven sys¬ 
tem yield a spin lifetime = 102 /is, more than twice 
the value measured for the {|0), | — l)}-qubit (40 fis). 

To examine the matter further, we inspect the ground 
state Hamiltonian of the NV centre, 

Hgs = 


/ 0 \ 

I Dg T d||F/2; T d_\_(^Ex T 1 

\ g^BpB d±{Ex - iEy) Dg + d\\Ez - Q^bB^ ) 

( 2 ) 

Here Dg ^ 2.87 GHz is the ground-state splitting, 
while d±_ = 0.17Hz m/V and d\\ = 35 x 10“^ Hz m/V 
are the perpendicular and parallel electric dipole com¬ 
ponents, respectively. = {iBx — By)/\/2 collects the 
transverse components of the magnetic field vector, with 
/3+ its complex conjugate. Bj and Ej^ with j = 
are magnetic fields or strain applied in the three spatial 
directions. We now examine the effect of these compo¬ 
nents individually. Transverse magnetic fields lead to 
an approximately quadratic dependence of the transition 
frequencies for Bx^y B^. However the dependence is 
^ 5 times smaller in the | ± 1) system compared to the 
|0) I ± 1) transitions at the bias field of 30 mT, since 
the shift caused by transverse magnetic fields is positive 
for both upper eigenstates, but negative for the lowest 
state. 


The effect of strain (or electric field fluctuations) E^ 
along the axis of the defect can be seen by inspection 
of the Hamiltonian to be a linear frequency shift for 
|0) ^ I 1) transitions, while causing no shift for the 
I — 1) 1 + 1) transition. The energies of the | ± 1) 

eigenstates at zero strain, Dg ± qiibBz are shifted to 

DgE^{giiBBzY + + ^y) under transverse strain. 

Such perturbations will therefore cause a frequency shift 
with twice the magnitude for the | ± 1) system. How¬ 
ever, due to the addition in quadrature, their effect is 
strongly suppressed in both manifolds by application of 
a bias field Bz. 

The improvement was observed for all NV centres ex¬ 
amined in the isotope-purified portion of the sample, but 
no improvement was found for NV centres in the portion 
with natural isotopic abundance. This distinction leads 
us to conclude that the that the longer coherence times 
measured for the | ± 1) may be brought about not by 
a better isolation from bath spins, but by a weaker de¬ 
pendence on environmental perturbations other than the 
measured magnetic field along the NV axis. 

We now describe how to make use of the improved co¬ 
herence lifetime in the | ± 1) system to measure arbitrary 
magnetic fields. The measurement principle is related to 
that used in the determination of the dephasing rate, as 
well as in previous magnetometry measurements using 
the {|0), |1)} and | ± 1) systems [17]. As a basis for com¬ 
parison, we measure the magnetic field sensitivity for the 
|0) ^ I — 1) transition using the pulsed electron spin reso¬ 
nance (ESR) method [Fig. [^a)] |T3]. The sensitivity can 
be found by recording the ESR spectrum of the transition 
[Fig. I b)] and determining the largest gradient Gmax-, 
so that 
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Gfnax 
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where tdc is the duration of the measurement sequence 
and N is the number of photons collected on average per 
sequence at the microwave detuning corresponding to the 
maximum gradient. 

We compare the sensitivity of the pulsed ESR measure¬ 
ment to a similar measurement in which the spin is coher¬ 
ently transferred along | — 1) ^ | + 1). In order to make 
optimal uses of the extended coherence time of the | ± 1) 
qubit, we transfer the system to | — 1) immediately after 
initialisation. The transition amplitudes for the effective 
two-level system are sensitive to the relative detuning 
Atot = A+ — A_ — 2giiB^Bz , with gfiB = 28 Hz/nT, and 
are therefore twice as sensitive as the simple qubit tran¬ 
sition. The longer coherence time of the {| + 1), | — 1)}- 
qubit allows for longer phase accumulation times, further 
increasing the sensitivity. 

We then measured the sensitivity for a more sophisti¬ 
cated sequence. Eor this measurement the spin is pre¬ 
pared in I — 1), then transferred into a superposition 
((a| +1) -h E'^I3\ — 1))/before being transferred 
further to the state | +1). In the experiment we use a nu¬ 
merically optimised pulse sequence provided by a genetic 
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FIG. 2. Sensing of arbitrary magnetic fields, a) Pulse sequence for pulsed ODMR in the {|0), | — l)}-manifold (bottom), a 
continuous two-tone | — 1) ^ | + 1) transfer pulse (middle) and a numerically optimised time-dependent pulse (top). The red, 
blue and green areas represent the Q+, and laser pulses applied. The 7r-pulses are performed on resonance, b) Detuning- 
dependent population for the sequences shown in a). c) Sensitivity as a function of sequence length for pulsed optically detected 
magnetic resonance (ODMR) (black), continuous transfer (orange) and optimised pulses (green), d) Measurement of the same 
arbitrary sub-/iT time-dependent magnetic held (dashed white line), compared for pulsed ODMR on the |0) ^ | — 1) transition 
(left), continuous transfer on the | — 1) ^ | + 1) transition (middle), and numerically optimised transfer on the | — 1) ^ | + 1) 
transition (right). The colour scale indicates the probability density for the magnetic held value, extracted from the NV 
huorescence for one second of integration time. 


algorithm [20] . The simulation takes into account the ro¬ 
tation of the state throughout the sequence to maximise 
the dependence on the magnetic held of the population in 
I — 1) after the sequence. The results show a reduction in 
the minimum detectable held, as shown in Fig. |c). Fi- 
nally, we show the effect of this improvement on the mea¬ 
surement of an actual varying magnetic held. The held 
is applied using a wire carrying a time-varying current /, 
which is spanned perpendicularly to the NV centre axis, 
thereby creating a magnetic held along the NV axis as de¬ 
picted in Fig. [^a). The held at each current value is mea¬ 
sured for one second, allowing us to extract a statistical 
uncertainty from the distribution of the recorded huores¬ 
cence counts. The greater sensitivity markedly improves 
this practical measurement, as shown in Fig. id). The 
techniques described above make possible the application 
of rephasing pulse sequences known from nuclear mag¬ 
netic resonance experiments and previous experiments 
with solid-state defects m- To measure the rephased co¬ 
herence time T 2 , we perform a similar sequence as for the 
T 2 measurement described above. This time however, 
after half the time delay, we rotate the state about the 
mean Bloch vector to perform a Carr-Purcell-Meiboom- 
Gill (CPMG) sequence. This rotation is achieved using 
synchronous pulses with Gaussian envelopes. Once 
again, using the | ± 1) states, the coherence time of the 


NV centre is increased to T 2 (±l) = (2.36 ± 0.09)ms , 
compared to T2(0,—1) = (1.66 dz 0.16) ms [Fig. |^a)]. 
While this increase is not as dramatic as for the T 2 values, 
it nonetheless represents a significant improvement, en¬ 
abling improved measurements of oscillating fields: The 
long phase accumulation time warranted by the rephas¬ 
ing sequence can be used to perform highly sensitive mag¬ 
netic field measurements, albeit only for known phases 
and frequencies of the magnetic field amplitude [21] . The 
smallest oscillating field that can be measured using this 
method is given by 


Ntac ^ ‘^Coexp[-{TAc/T2)P]SsgfiB ' 

Here N is the number of photons collected in a single 
acquisition and tac is the oscillation period of the field. 
Co is the normalized fluorescence contrast for the two 
states and p describes the shape of the coherence decay. 
Our measurements yield p = (2.03 ±0.38) for the {|0), | — 
1)} system andp = (2.02±0.1) for the |±1) system. The 
difference in magnetic quantum number Ss = 1 for the 
{|0), I — 1)} and 6s = 2 for the | ± 1) system. 

Once again, the factor-of-two increase in magnetic field 
dependence and the improved coherence time markedly 
reduce the minimum detectable field (see Fig. |b)). 
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FIG. 3. a) Comparison of T 2 coherence time measurements for the | zb l)-qubit (red, T 2 = (2.36 zb 0.09) ms) and the (|0), | — 1)) 
qubit (grey, T2(0, —1) = (1.66 zb 0.16) ms), b) Sensitivity to oscillating magnetic fields for both systems versus field oscillation 
period. The shaded areas indicate the sensitivity expected from the T 2 coherence time measurements, broadened by their 
uncertainty. Inset: Sinusoidal dependence of the fluorescence on the magnitude of an applied AC magnetic field for the two 
systems, recorded for phase accumulation times T 2 / 2 . c) Allan deviation for three measurement types: DC magnetometry with 
pulsed ODMR in the (|0), | — 1)) basis (black); DC magnetometry with optimized pulses in the (| — 1), | +1)) basis (green); and 
AC magnetometry in the (|0), | — 1)) basis (red). The dashed lines show the l/\/t dependence expected for purely shot-noise 
limited measurements. 


Measurements in the {|0), | zb l)}-system give a sensi¬ 
tivity of 6nT/Hz^/^, while significantly smaller changes 
in the field modulation amplitude of 2.2nT/Hz^/^ were 
detected using the | ± 1) system. The improvement is 
even more significant for longer field oscillation periods, 
as shown in Fig. [^b). 

A further important feature of the measurement tech¬ 
niques demonstrated here is the long-term stability given 
by the common-mode response of the | ± 1) levels to ex¬ 
traneous noise sources. Fig. c) shows the Allan devi¬ 
ation for DC and AC measurements. The measurement 
of static fields in the {|0), | — l)}-basis reaches its floor 
of 10 nT for time intervals of approximately one minute. 
This timescale is consistent with temperature fluctua¬ 
tions observed in the laboratory, which affect both the 
NV centre directly, as well as the permanent magnet used 
to apply the bias field. The same measurement performed 
using the optimised pulses for the |±1) qubit remains sta¬ 
ble for intervals which are an order of magnitude longer, 
and reaches a minimum deviation of ^ 800 pT. We pre¬ 
sume this sensitivity to be limited by actual fluctuations 
of the axial magnetic field around the sample. Neither of 
these mechanisms have a discernible effect on AC magne¬ 
tometry using the | ± 1) qubit, which displays shot-noise 
limited behaviour over several hours. The Allan devia¬ 
tion of this measurement reaches ^ 50 pT for intervals of 
approximately half an hour, limited only by the statistics 
and duration of our measurement. 

To conclude, the use of the full spin-1 manifold of the 


NV ground state gives access to numerous possibilities for 
storing and manipulating quantum states. The increase 
in coherence time given by encoding qubit states into the 
states I ± 1) when using isotopically purified diamond 
will be beneficial for quantum information applications, 
and leads to significant improvements in the sensitivity 
to static and modulated magnetic fields. The sensitivity 
is limited by the shot noise on the detected photons, and 
can be improved further by implementing existing micro- 
optical elements [22], or by using sparse NV ensembles. 
Potential improvements also include spin to charge con¬ 
version or photoelectron detection [23l [24] . 

We furthermore note that the use of two-photon ma¬ 
nipulation will enable the use of crossed-wire grids for 
the manipulation of arrays of NV centres, leading to 
a quadratic reduction in the number of connections re¬ 
quired for full control of qubit grids in large-scale quan¬ 
tum information systems [la usHm. Given these ad¬ 
vantages, we expect the techniques demonstrated in this 
work to be of broad use for applications of solid-state 
defects as sensors and carriers of quantum information. 
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